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We explore the influence of open-charm channels on charmonium properties, and profile the 1 3 D2, 
1 3 D3, and 2 x Pi charmonium candidates for X(3872). The favored candidates, the 1 3 D2 and 1 3 D3 
levels, both have prominent radiative decays. The 1 3 D2 might be visible in the D°D*° channel, 
while the dominant decay of the 1 3 D3 state should be into DD. We propose that additional discrete 
charmonium levels can be discovered as narrow resonances of charmed and anticharmed mesons. 
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I. INTRODUCTION 

Encouraged by the Belle Collaboration's sighting [1] 
of ?7^.(2 1 So) in exclusive B — » KKsK~tt + decays, we 
sketched a coherent strategy to explore rj' c and the re- 
maining charmonium states that do not decay into open 
charm, h c {l 1 I'i), 77 c2 (l 1 D 2 ), and ^ 2 (1 3 D 2 ), through B- 
meson gateways [2] . We argued that radiative transitions 
among charmonium levels and 7T7t cascades to lower-lying 
charmonia would enable the identification of these states. 

Now the Belle Collaboration has presented evidence [3] 
for a new narrow state, AT(3872) — > 7r + 7r~ Jfip, seen 
in B ± -> K ± X(3872). The CDF Collaboration has 
confirmed the new state in inclusive 1.96-TeV pp — ► 
A" (3872) + anything [4], as has the D0 experiment [5]. 
In addition, the CLEO [6], BaBar [7], and Belle [8] ex- 
periments have confirmed and refined the discovery of rj' c , 
fixing its mass and width as M(t]' c ) = 3637.7 ± 4.4 MeV 
and r{rf c ) = 19 ± 10 MeV [9]. 

In this Article, we develop the hypothesis that AT (3872) 
is a charmonium level. The new meson's position at 
D°D*° threshold makes it imperative to take account 
of the coupling between cc bound states and open- 
charm channels. Accordingly, we revisit the proper- 
ties of charmonium levels, using the Cornell coupled- 
channel model [10, 11] to assess departures from the 
single-channel potential- model expectations. Far below 
charm threshold, the nonrelativistic potential model is 
a good approximation to the dynamics of the charm- 
anticharm system. For excited states above the first few 
levels, the coupling of cc to charmed-meson pairs mod- 
ifies wave functions, masses, and transition rates. We 
estimate spin-splittings induced by communication with 
open-charm channels, and examine the effect of configu- 
ration mixing on radiative decay rates. We consider ip 2 



(1 3 D 2 ), ^3 (l 3 h>3), and h' c (2 1 P 1 ) as possible interpre- 
tations of AT (3872), commenting briefly on diagnostics 
of a general character that will help establish the na- 
ture of A" (3872). Independent of the identity of AT (3872), 
above-threshold charmonium states should be visible as 
narrow structures in 1 3 D 3 -> DD, 2 3 P 2 -» DD, DD*, 
1 3 F 4 -»• DD, DD*, and possibly 2 3 P -► DD. 

What do we know about A"(3872) ? Belle's clean sample 
of 36 events, entirely from B-meson decays, determines 
the mass of the new state as 3872.0 ± 0.6 ± 0.5 MeV, and 
yields a ratio of production x decay branching fractions, 



B(B+^K+X)B(X- 



0.063 ±0.014. (1) 



CDF observes 730 ± 90 events above background and 
determines a mass 3871.4 ± 0.7 ± 0.4 MeV. The ob- 
served mass lies 67 MeV above the l 3 Dj centroid in the 
potential- model template of Ref. [2]. The large number 
of events suggests that much of the CDF sample arises 
from prompt production of X(3872), not from B decays, 
and opens another path to the exploration of new char- 
monium states. Belle sets a 90% C.L. upper limit on the 
width, T(X(3872)) < 2.3 MeV. The tt+tt^ J/tp decay ap- 
pears to favor high dipion masses, but there is no detailed 
information yet about the quantum numbers J PC . Belle 
has searched in vain for radiative transitions to the l 3 Pi 
level; their 90% C.L. upper bound, 



T(V(3872) -> 7Xc i) 



r(AT(3872) 



< 0.89 



(2) 



conflicts with our single-channel potential-model expec- 
tations for the 1 3 D 2 state [2]. The theoretical estimate 
of the rnrJ/ip rate is highly uncertain, however. 



II. INFLUENCE OF OPEN-CHARM STATES 
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The Cornell group showed long ago that a very simple 
model that couples charmonium to charmed-meson de- 
cay channels confirms the adequacy of the single-channel 
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TABLE I: Thresholds for decay into open charm. 



Channel 



Threshold Energy (MeV) 



D°D° 
D+D- 
D°D*° or D*°D° 
D ± D *T 

D+D- 

D*°D*° 
D *+ D *- 

D+D*- or D*+D~ 

d:+d*- 



3729.4 
3738.8 
3871.5 
3879.5 
3936.2 
4013.6 
4020.2 
4080.0 
4223.8 



cc analysis below threshold and gives a qualitative un- 
derstanding of the structures observed above thresh- 
old [10, 11]. We now employ the Cornell coupled-channel 
formalism to analyze the properties of charmonium lev- 
els that populate the threshold region between 2M(D) 
and 2M (£>*), for which the main landmarks are shown 
in Table I. 

Our command of quantum chromodynamics is inade- 
quate to derive a realistic description of the interactions 
that communicate between the cc and cq + cq sectors. 
The Cornell formalism generalizes the cc model with- 
out introducing new parameters, writing the interaction 
Hamiltonian in second-quantized form as 

Ki = f ELi / : Pa(r)V(r r> (r') : d 3 r dV , (3) 

where V is the charmonium potential and p a ( r ) = 
^ip^ (r)X a ip(r) is the color current density, with V the 
quark field operator and A a the octet of SU(3) matrices. 
To generate the relevant interactions, ip is expanded in 
creation and annihilation operators (for charm, up, down, 
and strange quarks), but transitions from two mesons to 
three mesons and all transitions that violate the Zweig 
rule are omitted. It is a good approximation to neglect 
all effects of the Coulomb piece of the potential in (3). 

A full outline of the calculational procedure appears in 
Refs. [10, 11], but it is apt to cite a few elements here. 
We evaluate Eq. 3 between nonrelativistic (cc) states with 
wave functions determined by the Cornell potential, and 
l^o and l 3 Si cu, cd, and cs ground states with Gaussian 
wave functions. States with orbital angular momentum 
L > can decay in partial waves I = L =F 1. 

Following [10], we define a coupling matrix within the 
(cc) sector 



n nm (w) = 



E 



(n\n I \D l D J ){D l D ] \H I \m) 
(W - E Dt - E 3i + is) 



(4) 



where the summation runs over momentum, spin, and fla- 
vor. Above threshold (for W > Mi + Mj), £1 is complex. 
We decompose fl nm into a dynamical part (see [10]) that 
depends on the radial and orbital quantum numbers of 
the charmonium states and on the masses of Dj and Dj 
times the product rccoupling matrix shown in Table II 
that expresses the spin dependence for each partial wave. 



TABLE II: Statistical recoupling coefficients C, defined by 
Eq. D19 of Ref. [10], that enter the calculation of charmonium 
decays to pairs of charmed mesons. Paired entries correspond 
to l = L- \ and£ = L + l. 



State 


DD 


DD* 


D* D* 


'So 









: 2 




: 2 


3 Si 




i 




4 




7 




3 




3 




3 


3 Po 


1 : 





: 


: 


i 

3 


. 8 
' 3 


3 Pi 


U : 


U 


4 . 
3 ' 


. 2 
3 


U 


o 
: z 


x Pi 


: 





2 , 

3 


4 
3 


2 
3 


, 4 
J 3 


3 P 2 


: 


2 
5 


: 


6 
5 


4 . 
3 ' 


16 
15 


3 Di 


2 
3 


: 


2 
3 


: 


4 
15 


. 12 
' 5 


3 D 2 


: 





(i . 

5 


4 
' 5 


2 
5 


. 8 
J 5 


'Da 


: 





4 

5 ' 


6 
5 


4 
5 


6 

J 5 


3 D 3 


: 


3 
7 


: 


8 
7 


8 . 
5 • 


29 
35 


3 F 2 


3 
5 


: 


4 
5 


: 


11 

35 


, 16 
J 7 


3 F 3 


: 





8 . 
7 1 


fi 

' 7 


4 . 
7 ' 


10 
7 




: 





6 . 

7 ' 


8 

' 7 


6 . 

7 ' 


. 8 
J 7 


3 F 4 


: 


4 
9 


: 


10 

9 


12 
7 


. 46 
' 63 


3 G 3 


4 
7 


: 


6 
7 


: 


22 
63 


. 20 
' 9 


3 G 4 


: 





10 
9 


. 8 
' 9 


2 
3 


. 4 
' 3 


X G 4 


: 





8 . 

9 ' 


10 

9 


8 . 

9 ' 


10 
9 


3 G 5 


: 


5 
11 


: 


12 
11 


16 

9 


. 67 
' 99 



In each channel 2S+1 Lj, the physical states correspond 
to the eigenvalues of 



(5) 



The real parts of the energy eigenvalues are the char- 
monium masses. Imaginary parts determine the widths 
of resonances above threshold. The eigenvalues also de- 
termine the mixing among (cc) states and the overall 
fraction in the (cc) sector. 

To fix the (Coulomb + linear) charmonium potential, 



V(r) = —n/r + r/a 2 



(6) 



we adjust the strength of the linear term to reproduce 
the observed ip'-ip splitting, after including all the effects 
of coupling to virtual decay channels. Neglecting the in- 
fluence of open charm gives a = 2.34 GcV, k = 0.52, 
and a charmed-quark mass m c = 1.84 GcV. In the Cor- 
nell coupled-channel model, the virtual decay channels 
reduce the ip'-ip splitting by about 115 MeV, so the slope 
parameter has to be reduced to a = 1.97 GcV. 

The basic coupled-channel interaction (3) is spin- 
independent, but the hyperfine splittings of D and D* , 
D s and D*, induce spin-dependent forces that affect the 
charmonium states. These spin-dependent forces give 
rise to S-D mixing that contributes to the 0(3770) elec- 
tronic width, for example, and are a source of additional 
spin splitting, shown in the rightmost column of Ta- 
ble III. To compute the induced splittings, we adjust the 
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TABLE III: Charmonium spectrum, including the influence 
of open-charm channels. All masses are in MeV. The penul- 
timate column holds an estimate of the spin splitting due 
to tensor and spin-orbit forces in a single-channel potential 
model. The last column gives the spin splitting induced by 
communication with open-charm states, for an initially un- 
split multiplet. 



State 


Mass 


1% 


2 979.9° 


l 3 Si 


3 096.9° 


l 3 Po 


3 415.3° 


l 3 Pi 


3 510.5° 


l L Pi 


3 525.3 


1 3 P 2 


3 556.2° 


2 1 S 


3 637.7° 


2 3 Si 


3 686.0° 


l 3 Di 


3 769.9°' 


1 3 D 2 


3 830.6 




3 838.0 


1 3 D 3 


3 868.3 


2 3 P 


3 931.9 


2 3 Pi 


4007.5 


2 1 P 1 


3 968.0 


2 3 P 2 


3 966.5 



Centroid 



Splitting 
(Potential) 



Splitting 
(Induced) 



3 067.6" 



3 525.3 C 



3 673.9" 



(3 815)° 



3 968 d 



-90.5 
+30.2 

-114.9 C 
-11.6 C 
+1.5 e 
-31.9 C 

-50.4 
+16.8 

-40 




+20 

-90 
-8 

+25 



+2.8 
-0.9 

+5.9 
-2.0 
+0.5 
-0.3 

+15.7 
-5.2 

-39.9 
-2.7 
+4.2 
+19.0 

+10 
+28.4 
-11.9 
-33.1 



"Observed mass, from Review of Particle Physics, Ref. [13]. 



Inputs to potential determination. 
c Observed l 3 Pj centroid. 
d Computed. 

^Required to reproduce observed masses. 



bare centroid of the spin-triplet states so that the phys- 
ical centroid, after inclusion of coupled-channel effects, 
matches the value in the middle column of Table III. As 
expected, the shifts induced in the low- lying IS and IP 
levels are small. For the other known states in the 2S 
and ID families, coupled-channel effects are noticeable 
and interesting. 

In a simple potential picture, the ?7^(2 1 So) level lies 
below the ?//(2 3 Si) by the hyperfmc splitting given by 
M(ip') - M(rj' c ) = 327ra s |*(0)| 2 /9m2. Normalizing to 
the observed IS hyperfine splitting, M(J/ip) — M(t] c ) — 
117 MeV, we would find 



M(V>') - M(tQ = 67 MeV 



(7) 



which is larger than the observed 48.3 ± 4.4 MeV, as is 
typical for potential-model calculations. The 2S induced 
shifts in Tabic III draw ip' and rf c closer by 20.9 MeV, 
substantially improving the agreement between theory 
and experiment. It is tempting to conclude that the ip'-rj' c 
splitting reflects the influence of virtual decay channels. 

We lack a comprehensive theory of spin splittings for 
L > states, and various potential-model schemes differ 
appreciably in their predictions. (See Table I of Ref. [12] 
for a variety of estimates.) For the IP states, the spin 
splittings shown under Splitting (Potential) in Table III 



are those required to reproduce the observed masses; they 
are not predictions. For the ID and 2P levels, we have 
adopted as representative the spin splittings shown. 

To reproduce the observed mass of the l 3 Di ^(3770), 
we shift the bare ID centroid upward by 67.5 MeV. The 
other ID masses are thus pegged to the observed ^>(3770). 
In our model calculation, the coupling to open-charm 
channels increases the l 3 D2-l 3 Di splitting by about 
20 MeV, but does not fully account for the observed 
102 MeV separation between A(3872) and V(3770). It 
is noteworthy that the position of the 3 1 3 D 3 level 
turns out to be very close to 3872 MeV. For the 2P lev- 
els, we have no experimental anchor, so we adjust the 
bare centroid so that the 2 : Pi level lies at the centroid 
of the potential-model calculation. It is likely that we 
have more to learn about the influence of open-charm 
channels. 

The 2^1 level has been suggested [14] as an alterna- 
tive assignment for AT (3872) because it has an allowed 
7T7T transition to J/ip and a hindered Ml radiative tran- 
sition to the IP levels. The coupled-channel calculation 
places this state nearly 100 MeV above DD* threshold. 
As we shall see in quantitative detail presently, its al- 
lowed s-wave decay to D°D*° leads to an unacceptably 
large width, unless A (3872) lies below D°D*° threshold. 

The wave functions that correspond to physical states 
are linear combinations of potential-model cc eigenstates 
plus admixtures of charmed-meson pairs. We record the 
charmonium content of states of interest in Table IV. The 
open-charm pieces have the spatial structure of bound 
states of charmed mesons, but they are not molecular 
charm states in the usual sense: they are virtual con- 
tributions for states below threshold, and — unlike "deu- 
sons," for example [15] — they are not bound by one-pion 
exchange. 

Expectations for radiative transitions. As Table IV 
shows, the physical charmonium states are not pure 
potential-model eigenstates. To compute the El radia- 
tive transition rates, we must take into account both the 
standard (cc) — ► (cc)7 transitions and the transitions be- 
tween (virtual) decay channels in the initial and final 
states. Details of the calculational procedure are given 
in §IV.B of Ref. [11]. 

Our expectations for El transition rates among spin- 
triplet levels are shown in Table V. There we show both 
the rates calculated between single-channel potential- 
model eigenstates (in italics) and the rates that result 
from the Cornell coupled-channel model, to indicate the 
influence of the open-charm channels. The model repro- 
duces the trends of transitions to and from the Xc states 
in broad outline. Not surprisingly, the single-channel val- 
ues roughly track those calculated by Barnes & Godfrey 
in their potential [12]. For these low- lying states, the 
mixing through open-charm channels results in a mild 
reduction of the rates. 

We show the l 3 Di transition rates at the mass of 
V>(3770) and at the predicted l 3 Di centroid, 3815 MeV. 
For the V(3770), with its total width of about 24 MeV, 



4 



TABLE IV: Charmonium content of states near flavor threshold. The wave function $ takes account of mixing induced 
through open charm-anticharm channels. Unmixed potential- model eigenstates are denoted by \n 2s+1 Lj). The coefficient of 
the dominant eigenstate is chosen real and positive. The IS, IP, 2S, and 1 Di states are evaluated at their physical masses. 
The remaining ID and 2P states are considered at the potential-model centroids or at the masses in Table III?. We also 
show the 1 3 D2, 1 3 D3, and 2 x Pi states at the mass of X(3872). Z C c represents the (cc) probability fraction of each state. 

tf^So) = 0.986 |1 1 S > - 0.042 ^So) - 0.008 |3 1 S > - 0.002 ^So) - 0.001 |5 1 S >; Z cB = 0.974 
*(l 3 Si) = 0.983 |l 3 Si) - 0.050 |2 3 Si) - 0.009 |3 3 Si) - 0.003 |4 3 Si) + -0.001 |5 3 Si); Z cB = 0.968 

*(1 3 P ) = 0.919 |1 3 P ) - 0.067 |2 3 P ) - 0.014 |3 3 P ) - 0.005 |4 3 P ) - 0.002 j5 3 P ); Z cS = 0.850 
*(l 3 Pi) = 0.914 |l 3 Pi) - 0.075 |2 3 P!) - 0.015 |3 3 Pi) - 0.005 |4 3 P!) - 0.002 |5 3 Pi); Z cB = 0.841 
tf^Pi) = 0.918 I^Pi) - 0.077 |2 1 P 1 ) - 0.015 ft 1 Pi) - 0.005 |4 1 Pi> - 0.002 js^i); Z c - C = 0.845 
*(1 3 P 2 ) = 0.920 |1 3 P 2 ) - 0.080 |2 3 P 2 ) - 0.015 |3 3 P 2 ) - 0.005 j4 3 P 2 ) - 0.002 |5 3 P 2 ) - 0.002 |1 3 F 2 ); Z cB = 0.854 

*(2 1 S ) = 0.087 |l 1 So> + 0.883 ^S,,) - 0.060 |3 1 S > - 0.016 |4 1 S ) - 0.007 |5 1 S ) - 0.003 |6 1 S >; Z a - C = 0.791 
*(2 3 Si) = 0.103 |l 3 Si) + 0.838 |2 3 Si) - 0.085 |3 3 Si) - 0.017 |4 3 Si) - 0.007 |5 3 Si) - 0.002 j6 3 D!) 
+0.040 |l 3 Di) - 0.008 |2 3 Di); £ cg = 0.723 

*(l 3 Di) = 0.694 | l 3 Di) +0.097 e ' 9351 " |2 3 Di) + 0.008 e" ' 668 ^ |3 3 Di) + 0.006 e ' 904 " 1 |4 3 D!) 

+0.013 e 7 * 2 ™ |l 3 Si) + 0.168 e 805 ^ |2 3 Si) + 0.014 e - 866 ^ |3 3 Si) + 0.012 e -°- 22 «* |4 3 Si) 
+0.001 e °- 27ai7V |5 3 Si> + 0.001 e -°- 267in |6 3 Si); Z cB = 0.520 
*(1 3 D 2 ) = 0.754 |1 3 D 2 ) - 0.084 |2 3 D 2 ) - 0.011 |3 3 D 2 ) - 0.006 |4 3 D 2 >; Z c£ = 0.576 
*(1 1 D 2 ) = 0.770 |1 1 D 2 ) - 0.083 |2 1 D 2 ) - 0.012 |3 1 D 2 ) - 0.006 |4 1 D 2 ); Z cS = 0.600 
*(1 3 D 3 ) = 0.812 1 1 3 D 3 ) +0.086 e " ^ |2 3 D 3 ) + 0.013 e" ' 969 ^ |3 3 D 3 ) + 0.007 e ' 980 " 1 |4 3 D 3 ) 

+0.016 e ' 848 ^ |1 3 G 3 ) + 0.003 e- ' 291 " 1 |2 3 G 3 ); Z cB = 0.667 

*(2 3 P ) = 0.040 e" 454 ^ |1 3 P > +0.532 |2 3 P > + 0.024 e -°- 889i7r |3 3 P ) + 0.010 e ' 867 " 1 |4 3 P ) 

+0.006 e -°- 976i7r |5 3 P ); Z cB = 0.286 
*(2 3 P!) = 0.218 e -° A56in |l 3 Pi) +0.821 |2 3 Pi) +0.058 e ' 516 ^ |3 3 Pi> + 0.032 e °- 976i * |4 3 Pi) 

+0.008 e °- 986i7r |5 3 Pi); Z cS = 0.726 
*(2 1 Pi) = 0.216 e" ' 226 ^ I^Pi) +0.852 |2 4 Pi) + 0.079 e°- 780<,r ^Pi) + 0.023 e" 890 ^ |4 x Pi) 

0.007 e 985 ™ IS 1 ?!); Z ce = 0.883 
4>(2 3 P 2 ) = 0.234 e" 046 ^ |1 3 P 2 ) + 0.754 |2 3 P 2 ) + 0.097 e ' 876 ^ |3 3 P 2 ) + 0.016 e" ' 743 ^ |4 3 P 2 ) 

0.007 e °- 898in |5 3 P 2 ) +0.370 e ' 775 " 1 |1 3 F 2 ) + 0.035 e -°- 317 ™ |2 3 F 2 ) + 0.002 e 097 ™ |3 3 F 2 ); Z c - C = 0.771 

M = 3872 MeV : *(1 3 D 2 ) = 0.596 |1 3 D 2 ) - 0.108 |2 3 D 2 ) - 0.004 |3 3 D 2 ) - 0.006 |4 3 D 2 ); Z cB = 0.367 

M = 3872 MeV : *(1 3 D 3 ) = 0.813 |1 3 D 3 ) + 0.089 e ' 989 ^ |2 3 D 3 ) + 0.013 e" 965 ^ |3 3 D 3 ) + 0.007 e °- 978in |4 3 D 3 > 

+0.017 e ' 837 " 1 |1 3 G 3 ) + 0.003 e^ ' 305 " 1 |2 3 G 3 ); Z cS = 0.669 
M = 3872 MeV : *(2 1 Pi) = 0.134 e -0 004 ^ |l 1 Pi) + 0.374 ^Pi) + 0.035 e ' 993 ^ ^Pi) + 0.003 e" ' 981 ^ |4 4 Pi) 
+0.004 e ' 996 *" ^Pi); £ cg = 0.159 



the l 3 Di(3770) — > Xco 7(338) transition might someday 
be observable with a branching fraction of 1%. 

For the 1 3 D 2 and 1 3 D 3 levels, we have computed 
the radiative decay rates at the predicted l 3 Di cen- 
troid, 3815 MeV, at the mass calculated for the states 
(3831 MeV and 3868 MeV, respectively), and at the mass 
of A" (3872). We will compare the partial widths for the 
Xci 7(344) and x C 2 7(303) with the expected ir + n~J/ip 
and open-charm decay rates presently. 

We have evaluated the radiative decay rates for the 
2 3 Pj levels at the calculated centroid and at the pre- 
dicted mass, where that is displaced appreciably from 
the centroid. We shall see below that all of these rates 
are small compared to the expected open-charm decay 
rates. 

Expectations for hadronic transitions. The Beijing 
Spectrometer (BES) observation [16] of a branching frac- 
tion £(V>(3770) -> 7T+7T- J/ip) = (0.59 ± 0.26 ± 0.16)% 
would imply a hadronic cascade rate r(l 3 Di — ► nnJ/ip) « 
210 ± 130 keV, considerably larger than the 45 kcV, in- 



ferred [17] from older data, that we took as normaliza- 
tion in Rcf. [2]. By the Wigner-Eckart theorem for El-El 
transitions, all the l 3 Dj — > nnJ/ip rates should be equal 
(for degenerate 3 Dj states), so the higher BES normal- 
ization would increase r(l 3 D 2 — > irirJ/ip) to hundreds 
of keV, as remarked by Barnes & Godfrey [12]. Com- 
bined with our estimate that r(l 3 D 2 (3872) — > 7Xci) ~ 
207 keV, the larger irirJ/ip rate relaxes somewhat — but 
does not eliminate — the tension between the 3 D 2 assign- 
ment for X and the Belle bound in Eq. (2). 



The BES rate, which is based on a handful of events, 
is challenged by a CLEO-c limit [9], £(V>(3770) -> 
TT+n-J/tp) < 0.26% at 90% C.L. Both experiments are 
accumulating larger data samples that should improve 
our knowledge of this important normalization. 
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TABLE V: Calculated and observed rates for El radiative 
transitions among charmonium levels. Values in italics result 
if the influence of open-charm channels is not included. 



Transition 

(7 energy in MeV) 



Partial width (keV) 
Computed Measured 



Xco J/0 7(303) 
Xd^J/V- 7(390) 
Xc2^ J/4> 7(429) 

V/^Xc2 7(129) 
1>' - Xd 7(172) 
V/^XcO 7(261) 



3770) — X c2 7(208) 
3770) -> Xci 7(251) 
3770) ^Xco 7(338) 

3815) ^ X c2 7(250) 
3815) -»xd 7(293) 
3815) -» Xco 7(379) 

3815) ^Xc2 7(251) 
3815) -> Xci 7(293) 

3831) ^x C 2 7(266) 
3831) -* Xcl 7(308) 

3872) -> x C 2 7(303) 
3872) - Xc i 7(344) 

3815) ~> Xc 2 7(251) 
3868) -► Xc2 7(303) 
3872) x C 2 7(304) 

3933) -f J/i/> 7(747) 
3933) -> V' 7(239) 
3933) -f t/)(3770) 7(160) 

3968) -»■ J/i/) 7(775) 
3968) -f V' 7(272) 
3968) -»■ V(3770) 7(193) 

3968) -»■ J/0 7(775) 
3968) -f V' 7(272) 
3968) -»■ 0(3770) 7(193) 
3968) -f 1 3 D 2 (3815) 7(150) 
3968) -»■ 1 3 D 2 (3831) 7(135) 
3968) -f 1 3 D 2 (3872)7(95) 

4012) J/0 7(811) 
4012) -f 0' 7(313) 
4012) -f ^(3770) 7(235) 

3968) -> J/0 7(775) 
3968) -f V' 7(272) 
3968) -»■ 0(3770) 7(193) 
3968) -f 1 3 D 2 (3815) 7(150) 
3968) -»■ 1 3 D 2 (3835) 7(131) 
3968) 1 3 D 2 (3872)7(95) 
3968) -»■ 1 3 D 3 (3815) 7(150) 
3968) -f 1 3 D 3 (3868) 7(99) 
3968) 1 3 D 3 (3872)7(95) 



iJ5 -> 107 


119±19° 


22c? -> 216 


291±48 a 


300 -> 287 


426±51 a 


^3 -> 23 


27±4 6 


33 -» 32 


27±3 6 


36 -> 38 


27±3 6 


3.2 -»■ 3.9 




i«3 -f 59 




25^ -> 225 




5.5 -*■ 6.8 








3^ -> 371 




50 -> 40 




230 -» 191 




55 — > 45 




2&£ -> 212 




— > 4o 




362 -> 207 








325 286 




3^ J -> 299 




95 — ► 19 




1Z 1 — > 6<$ 




$g _ > 11 




i.Z0 — ► 77 




ic?0 — ► 155 




iOJ -> 43 




iiO — ► 68 




JSC — ► 102 




25 — > 5 




37 — > 1.8 




25 — > 0.25 




i0 — > 0.23 




i32 ->■ 94 




260 -> 151 




^3 -> 11 




-> 19 




ic?0 -> 314 




i.O -»■ 1.4 




7.^ -> 18 




5 -» 12 




i.P 3.4 




-> 82 




i2 26 




-> 23 





"Derived from the 2003 "unchecked fit" of Ref. [13]. 
'Branching fractions from CLEO via Skwarnicki [9]. 



III. DECAYS INTO OPEN CHARM 

The calculated partial widths for decays of charmo- 
nium states into open charm appear in Table VI. Ex- 
perience [10] teaches that once the position of a res- 
onance is given, the coupled-channel formalism yields 
reasonable predictions for the other resonance proper- 
ties. The l 3 Di state V"(3770), which lies some 40 MeV 
above charm threshold, offers an important benchmark: 
we compute r(^"(3770) -» DD) = 20.1 MeV, to be 
compared with the Particle Data Group's fitted value 
of 23.6 ± 2.7 MeV [13]. The variation of the l 3 Di width 
with mass is shown in the top left panel of Figure 1. 

Barnes & Godfrey [12] have estimated the decays of 
several of the charmonium states into open charm, using 
the 3 Po model of qq production first applied above charm 
threshold by the Orsay group [18]. They did not carry 
out a coupled-channel analysis, and so did not determine 
the composition of the physical states, but their estimates 
of open-charm decay rates can be read against ours as a 
rough assessment of model dependence. 

The long-standing expectation that the 1 3 D 2 and 1 X D 2 
levels would be narrow followed from the presumption 
that these unnatural parity states should lie between 
the DD and DD* thresholds, and could not decay into 
open charm. At 3872 MeV, both states can decay into 
D°D*°, but the partial widths (Table VI) are quite small. 
We show the variation of the 1 3 D 2 partial width with 
mass in the top right panel of Figure 1; over the re- 
gion of interest, it does not threaten the Belle bound, 
T(A(3872)) < 2.3 MeV. The range of values is quite 
similar to the range estimated for r(l 3 D 2 — ► mrJ/ip), so 
we expect roughly comparable branching fractions for de- 
cays into D°D*° and tt+tt" J/ip. If X(3872) docs turn out 
to be the 1 3 D 2 level, we expect M(1 1 D 2 ) = 3880 MeV 
and r(l 1 D 2 -> D°D*°) w 1.7 MeV. 

The natural-parity 1 3 D3 state can decay into DD, but 
its /-wave decay is suppressed by the centrifugal bar- 
rier factor, so the partial width is less than 1 MeV at a 
mass of 3872 MeV. Although estimates of the hadronic 
cascade transitions are uncertain, the numbers in hand 
lead us to expect r(l 3 D 3 -> tt+tt-J/^) < ±r(l 3 D 3 -> 
DD), whereas r(l 3 D 3 -> 7 x C 2) ~ |r(l 3 D 3 -> DD), if 
A (3872) is identified as 1 3 D 3 . The variation of T(1 3 D 3 -> 
DD) with mass is shown in the middle left panel of Fig- 
ure 1. Note that if 1 3 D 3 is not to be identified with 
A (3872), it may still be discovered as a narrow DD res- 
onance, up to a mass of about 4000 MeV. 

In their study of B+ -» AT+V>(3770) decays, the Belle 
Collaboration [19] has set 90% CL upper limits on the 
transition B+ -> A^+A(3872), followed by A (3872) -> 
DD. Their limits imply that 



S(A(3872) 
£(A(3872) - 



D°D°) 



< 48(A 
D+D~) < 3B(X 



-JH>) , 

-m ■ 



(8) 



This constraint is already intriguingly close to the level 
at which we would expect to see 1 3 D 3 — > DD. 
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The constraint on the total width of V(3872) raises 
more of a challenge for the candidate, whose s- 

wave decay to D°D*° rises dramatically from threshold, 
as shown in the middle right panel of Figure 1. Within 
the current uncertainty (3871.7 ± 0.6 MeV) in the mass 
of X, the issue cannot be settled, but the 2 : Pi interpre- 
tation is viable only if X lies below D°D*° threshold. If 
a light 2 1 P 1 does turn out to be AT(3872), then its 2 3 Pj 
partners should lie nearby. In that case, they should be 
visible as relatively narrow charm-anticharm resonances. 
At 3872 MeV, we estimate r(2 3 P! -> DD*) « 21 MeV 
and r(2 3 P 2 -> DD) w 3 MeV. The bottom left panel 
in Figure 1 shows that the 2 3 P2 level remains relatively 
narrow up to the opening of the D*D* threshold. 

The 2 3 P state is an interesting special case, as illus- 
trated in the bottom right panel of Figure 1. Through 
the interplay of nodes in the radial wave function, 
form-factor effects, and the opening of new channels, 
r(2 3 P -» DD) decreases from w 60 MeV at 3872 MeV 
to about 12 MeV near 3930 MeV. The total width for 
decay to open charm then rises in steps as M(2 3 P ) in- 
creases through the D S D S and D*D* thresholds. To es- 
timate the competing annihilation decay rate, we scale 
r(2 3 P — > gg — > hadrons) w r(l 3 P — > gg — > hadrons) • 
|i?^ P (0)| 2 / |i?' 1P (0)| 2 , where R'(0) is the derivative of the 
radial wave function at the origin. This yields r(2 3 Po — > 
gg -> hadrons) « 1.36 x 10.6 MeV = 14.4 MeV [13], an 
estimate that should probably be reduced by the |2 3 Po) 
fraction of the physical 2 3 Po state. 

We call attention to one more candidate for a narrow 
resonance of charmed mesons: The 1 3 F4 level remains 
narrow (r(l 3 F 4 charm) < 5 MeV) up to the _D*D* 
threshold. Its allowed decays into DD and DD* are 
inhibited by £ = 4 barrier factors, whereas the D*D* 
channel is reached by I = 2. 



IV. FOR THE FUTURE 

On the experimental front, the first order of busi- 
ness is to establish the nature of V(3872). Determin- 
ing the spin-parity of X will winnow the field of can- 
didates. The charmonium interpretation and its promi- 
nent rivals require that V(3872) be a neutral isoscalar. 
Are there charged partners? A search for V(3872) — > 
7T O 7r°J/0 will be highly informative. As Barnes & God- 
frey [12] have remarked, observing a significant 
signal establishes that X is odd under charge conju- 
gation. Voloshin has commented [20] that the ratio 
TZq = T(X -> 7r°7r°J/V')/r(V -> n+n-Jfip) measures the 
dipion isospin. Writing Ti = T(X — > (7r + 7r~)/ J/0), we 
see that TZo — ^/(l+ri/r ), up to kinematic corrections. 
Deviations from TZo — \ signal the isospin-violating de- 
cay of an isoscalar, or the isospin-conserving decay of 
an isovector. Radiative decay rates and the prompt (as 
opposed to -B-decay) production fraction will provide im- 
portant guidance. Other diagnostics of a general nature 
have been discussed in Refs. [12, 14, 21, 22]. 



TABLE VI: Partial widths for decays of charmonium states 
into open charm, computed in the Cornell coupled-channel 
model. All masses and widths are in MeV. Only significant 
partial widths are tabulated, and "total" refers to the sum 
of open-charm decays. Properties of the candidate states for 
X(3872) and their partners are evaluated at 3872 MeV and 
also at the potential-model centroid for each state. Decays 
occur with orbital angular momentum I. For DD* modes, 
the sum of DD* and DD* is always implied. 



State 


Mass 


£ 


Channel 


Width 


Total Width 


1 Di 


3770 


1 


D°D° 
D+D- 


11.8 
8.3 


20.1 


1 3 D 2 


3815 












1 3 D 2 


3872 


1 


D°D*° 


0.045 


0.045 


Pd 2 


3815 












Pd 2 


3872 


1 


D°D*° 


0.030 


0.030 


Pd 2 


3880 


1 


D°D*° 


1.7 


1.7 






Q 
O 


D°D° 
D+D- 


0.47 
0.39 


u.ou 


1 3 D 3 


3902 


3 


D°D° 
D+D- 


0.84 
0.72 


1.56 


2 3 P 


3872 





D"D" 
D+D- 


27 
32 


59 


2 3 P 


3930 





D°D° 
D+D' 

D°D° 


5.0 
7.4 

0.27 


12.4 


2 3 P 


3968 





D+D~ 
D S D S 


0.85 
40 


41.1 


2 3 Pi 


3872 





D°D*° 


20.9 


20.9 


2 3 Pi 


3968 





D°D*° 
D+D*- 


71.4 
78.9 


150.3 


2 x Pi 


3871.6 





D D 


4.28 


4.28 






u 


D°D*° 
D+D*- 


35.5 
39.2 


7/1 7 


^ 12 


OO ( z 


o 
Z 


D°D° 
D+D- 

D°D° 


1.63 
1.42 

4.4 


O.UO 


2 3 P 2 


3968 


2 


D+D- 
D° D*° 
D+D*- 


4.2 
2.57 
2.16 


13.4 








D°D° 


46 










D+D- 


45 




1 3 F 2 


4054 


2 


D S D S 
D°D*° 
D+D*- 


4 
31 
29 


155 


1 3 F 3 


4054 


2 


D°D*° 
D+D*- 


44.9 
42.6 


87.5 


Pf 3 


4054 


2 


D°D*° 
D+D*- 


32.2 
30.8 


66.4 








D°D° 


1.96 




1 3 F 4 


4054 


4 


D+D- 
D°D*° 
D+D*- 


1.78 
0.62 
0.50 


4.88 
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Within the charmonium framework, A (3872) is most 
naturally interpreted as the 1 3 D 2 or 1 3 D 3 level, both of 
which have allowed decays into irnJ/ip. The 2 1 3 D2 
state is forbidden by parity conservation to decay into 
DD but has a modest D°D*° partial width for masses 
near 3872 MeV. Although the uncertain mrJ/ip partial 
width makes it difficult to estimate relative branching 
ratios, the decay A(3872) — ► Xd 7(344) should show 
itself if X is indeed 1 3 D2. The \c2 7(303) line should 
be seen with about j the strength of \ci 7(344). In 
our coupled-channel calculation, the 1 3 D 2 mass is about 
41 MeV lower than the observed 3872 MeV. In contrast, 
the computed 1 3 D 3 mass is quite close to 3872 MeV, and 
1 3 D 3 does not have an El transition to Xci 7(344). The 
dominant decay of the 3~~ 1 3 D3 state should be into 
DD; a small branching fraction for the mrJ/ip discovery 
mode would imply a large production rate. One radia- 
tive transition should be observable, with r(A(3872) — > 
Xc2 7(303)) > r(A(3872) -> n + n-Jfi/>). We under- 
score the importance of searching for the \ci 7(344) and 
Xc2 7(303) lines. 

Beyond pinning down the character of A (3872), ex- 
periments can search for additional narrow charmonium 
states in radiative and hadronic transitions to lower- 
lying cc levels, as we emphasized in Ref. [2], and in neu- 
tral combinations of charmed mesons and anticharmcd 
mesons. The coupled-channel analysis presented in this 
paper sets up specific targets. 

On the theoretical front, we need a more complete un- 
derstanding of the production of the charmonium states 
in B decays and by direct hadronic production, includ- 
ing the influence of open-charm channels. Understand- 
ing of the production mechanisms for molecular charm 
or ccg hybrid states is much more primitive. We need to 
improve the theoretical understanding of hadronic cas- 
cades among charmonium states, including the influence 
of open-charm channels. The comparison of charmonium 



transitions with their upsilon counterparts should be in- 
formative. The analysis we have carried out can be ex- 
tended to the bb system, where it may be possible to see 
discrete threshold-region states in direct hadronic pro- 
duction. Because the Cornell coupled channel model is 
only an approximation to QCD, it would be highly desir- 
able to compare its predictions with those of a coupled- 
channel analysis of the 3 Po model of quark pair produc- 
tion. Ultimately, extending lattice QCD calculations into 
the flavor-threshold region should give a firmer basis for 
predictions. 

In addition to the UPi h c , the now-established 
2 1 So r]' c , and the long-sought 1 1 D 2 t] c2 and 1 3 D 2 V2 
states, discrete charmonium levels are to be found as nar- 
row charm- ant icharm structures in the flavor-threshold 
region. The most likely candidates correspond to the 
1 3 D 3 , 2 3 P 2 , and 1 3 F 4 levels. If X(3872) is indeed a char- 
monium state — the 3 D 2 and 3 D 3 assignments seem most 
promising — then identifying that state anchors the mass 
scale. If A (3872) is not charmonium, then all the char- 
monium levels remain to be discovered. Finding these 
states — and establishing their masses, widths, and pro- 
duction rates — will lead us into new terrain. 
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FIG. 1: Partial and total widths near threshold for decay of charmonium states into open charm, computed in the Cornell 
coupled-channel model. Long dashes: D°D°, dots: D + D~ , dot-dashes: D°D*°, dashes: D + D*~ , thin line: D*°D*°, short 
dashes: D* + D*~ , widely spaced dots: D S D S , thick line: sum of open-charm channels. Belle's 90% C.L. upper limit [3], 
T(X(3872)) < 2.3 MeV, is indicated on the window. For DD* modes, the sum of DD* and DD* is always implied. 



